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Abstract 



We consider the scattering of an electron from a semi-infinite one- 
dimensional random medium. The random medium is characterized by force, 
—dV/dL being the basic random variable. We obtain an analytical expression 



a 
o 

for the stationary delay time (r) distribution P s (t) within a random phase 



approximation. Our result agrees with earlier analytical expressions, where 

> 

0^ ' the random potential is taken to be of different kind, indicating universality 

^" ! 

of the delay time distribution, i.e., delay time distribution is independent of 

^ ' the nature of disorder. 
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In recent years universal parametric correlations of phase shifts and delay times in meso- 
scopic systems is being studied intensively |jlj]. The time delay in a scattering event (or 
duration of a collision event) is an interesting aspect in itself in the general theory of quan- 
tum scattering. Wigner was the first to establish the relation between the time delay and 
the energy derivative of scattering phase shift |2|. Distribution of delay times in quantum 
chaotic regime have been shown to be universal as it depends only on the symmetry prop- 
erty of the Hamiltonian or scattering matrix [|l]|§. The delay time statistics is intimately 
connected with the issue of dynamic admittance of microstructres (or mesoscopic systems) 
lUI]], for example quantum capacitance and its fluctuation ||. The wave packet incident on 
the surface a sample is not backscattered (or reflected) immediately. There will be some time 
delay before it is reflected. This leads to a non-cancellation of the instantaneous currents at 
the surface involving the incident and the reflected wave. This in-turn is expected to lead to 
a low temperature l//-type noise for the fluctuating surface currents in the random systems 
0-§|]. The study of change of density of states due to scatterer is also directly related to the 
phase derivative of scattering phase shift with respect to the energy, i.e., to the delay time. 

The distribution of delay time and its correlations in higher dimensions, where system 
exhibits the Anderson localization has not been addressed so far. The first study of the 
stationary distribution P s (t) of delay time r for a disordered semi-infinite one-dimensional 
chain was carried out in reference |7j]. Here authors used the invariant imbedding approach. 
The underlying random potential V(x) is treated as a Gaussian white noise with zero mean. 
Using the random phase approximation (RPA) analytical expression for the P s (t) was ob- 
tained, which exhibits 1/r 2 dependence for the tail of the delay distribution. Further de- 
velopments for P s (t) using supersymmetric potentials lead to same distribution function 
P s (t) within RPA. This has lead to a conjecture that within RPA, P s (t) is independent 
of nature of disorder and, hence, is universal |J. Our recent numerical study has clearly 
indicated that long time delay distribution is universal beyond RPA JTIJ. In our present 
work we calculate the distribution of delay time where we take —dV/dx as the basic random 
variable with delta correlated Gaussian distribution and we obtain analytical expression for 



P s {t) in RPA. The stationary distribution has the same functional form obtained earlier 
with different random potential indicating the universal nature of P s (t). 
The model Hamiltonian for the 1-D disordered system is 

-h 2 d 2 

where V(x) for < x < L is the random potential. The disordered sample extends from 
x = to X = L, the two ends being connected to perfect leads. Consider an electron of wave 
number k incident at x = L from right. It is partially reflected with the complex reflection 
coefficient R(L) and partially transmitted. The transmission and reflection coefficients are 
emergent quantities of direct physical interest for the conductance problem. The method 
of invariant imbedding was proposed originally by S. Chandrashekhar in the context of 
radiative transfer through stellar atmosphere [JO]]. His method consists of viewing the given 



sample of length L as imbedded in a larger sample of length L + AL and then setting up 
an equation for the resulting change in the S- matrix AS as AL — > 0. In order to look for 
the complex reflection coefficient, we transform Eqn. [I] to the invariant imbedding equation 
for the complex reflection amplitude R(L) = \R(L)\Exp(i6(L)) . The evolution equation for 
R{L) is now given by |T2j 

dR(L) 



dL 



with 



h{L) + 2if a (L)R(L) - h(L)R'(L), (2) 
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k{L)dL' 
fo(L) = k(L) 



and 
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k 2 = —{E-V{L)). 



The above equation Eqn. || was first studied in Ref. |H| to evaluate the resistance and 



its fluctuation in a disordered quantum wire. The invariant imbedding method has been 
generalized to N-channel case and in an equivalent random phase approximation has lead 



to DMPK (Dorokhov-Mello-Pereyra-Kumar) equation [14]], using which coherent transport 
properties have been analyzed extensively in mesoscopic systems. 

In the present problem we consider random potential V(L) to be bounded having a small 
amplitude. However, —dV/dL can be unbounded and we treat £,{L) = —dV/dL as our basic 
random variable. The energy of incident electron is assumed to be large, i.e., much larger 
than the magnitude of the upper bound on the potential. In that case we have 



/l(L) = If ^ and fo(L) = V f^ (3) 

We take £,{L) to be Gaussian delta correlated random number with zero mean and 

(am(L')) = 2a5(L-L'). (4) 

Here, the < .... > denotes the ensemble average with respect to the realizations of the 
stochastic variable £ and a denotes the strength of the disorder. The equation for the phase 
(6) is readily obtained from Eqn. as 



where we have set \R\ — 1 since we will be interested in the limit L — > oo (semi- 
infinite medium), i.e., total back-reflection with probability one. The delay time is given by 
t = HdO/dE. Differentiation of Eqn. |5] with respect to E leads to the following equation for 
the evolution of r: 

dr \j2m 2ft /r . . ,„. 2 ., rS , nS . . 

+ —S(L)sin(6) - -£(L)cos(9)t (6) 



dL ^E £ 2 E 

From Eqns. [5] and |^ one can obtain readily obtain the equation governing the evolution of 



the joint probability distribution W(r, 9; L) for 6 and r by using the Van Kampen lemma [IS 
and Novikov's theorem p~6|— 118|| . In our case, however, we are interested only in the marginal 
probability distribution P(t;L) = J n W(t,6; L)d9 of delay time r. The delay time being 
the derivative of phase we expect it to fluctuate much more rapidly as compared to the 
phase itself. We therefore make the decoupling approximation, as done in earlier literature 



by |7|||, treating 9 and r as statistically independent variables in the large length (L) limit. 
As mentioned earlier, we are interested in the case of high energy particles (E ^> V) and 
in this limit the distribution of 9 becomes uniform |19|-pT|, i.e., P(9) = 1/2-k. This is 
generally referred to as random phase approximation (RPA). Within the above mentioned 
approximations after a straight forward algebra the evolution equation for P(t; L) becomes 

dP „ d (2ahdP V2^ „ 4a „ 2a 9 dP) 

= ft — < P A tP -\ t > 7) 

8L dr\ E* dr ^E hE 2 TiE 2 8t J K J 

The stationary distribution P s {t) for r in the limit L — > oo can be obtained by setting 
dP/dL = 0. We get the following expression for normalized P s {t) 

\ p \tan~ 1 T 

«W = (, W - 1)|1+t!| ( 8 ) 

In the above expression we have redefined r in a dimensionless form r = rE/h and 
A = \JlmE E 2 /(2ah). The most probable value of r occurs at r = A/2. As r — > oo, 
P s (t) —> 1/t 2 , i.e., the distribution has a long time tail which goes as 1/t 2 . This leads 
to the logarithmic divergence of the average value of r indicating that the origin of such 
a tail is due, presumably, to the Azbel resonances E| which make Landauer's four probe 



conductance infinite even for a finite sample. In case of these resonant realizations,the 
time spent by the particle inside the medium is large as it travels a large distance before 
getting reflected. It is now well established that coherent interference effects, due to elastic 
scattering by the serial static disorder lead to localization of eigenstates for arbitrary weak 
disorder. The localization length I of these eigenstates is a self averaging quantity [23|] and 



in a one- dimensional system it is directly proportional to elastic mean free path. The most 
probable value r max of r is proportional to a time taken by a particle to traverse a distance 
of the order of localization length /, T max oc 2l/(hk/m), where k is the incident energy. From 
this one can readily obtain the behavior of localization length on the material parameters, 
namely, I ex E 2 /a. 

Our above analytical expression has the same functional form as obtained earlier where 
potential itself is treated as a Gaussian random variable using a different invariant imbedding 



equation. From this we conclude that two different models of random variable lead to the 
same universal distribution of the delay time. Thus reinforcing the conjecture of universal 
behavior of the delay time distribution, independent of nature of disorder within RPA. 
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